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Abstract 

In this paper, we have considered the mechanical stability of a jellium system 
in the presence of spin degrees of freedom and have generalized the stabilized 
jellium model, introduced by J. P. Perdew, H. Q. Tran, and E. D. Smith 
[ Phys. Rev. B 42, 11627 (1990)], to a spin-polarized case. By applying this 
generalization to metal clusters (Al, Ga, Li, Na, K, Cs), we gain additional 
insights about the odd-even alternations, seen in their ionization potentials. 
In this generalization, in addition to the electronic degrees of freedom, we 
allow the positive jellium background to expand as the clusters' polarization 
increases. In fact, our self-consistent calculations of the energetics of alkali 
metal clusters with spherical geometries, in the context of density functional 
theory and local spin density approximation, show that the energy of a cluster 
is minimized for a configuration with maximum spin compensation (MSC). 
That is, for clusters with even number of electrons, the energy minimization 
gives rise to complete compensation (JVf = N\), and for clusters with odd 
number of electrons, only one electron remains uncompensated (N^ — N^ = 1). 
It is this MSC-rule which gives rise to alternations in the ionization potentials. 
Aside from very few exceptions, the MSC-rule is also at work for other metal 
culsters (Al, Ga) of various sizes. 
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I. INTRODUCTION 



The subject of metal clusters has gained considerable momentum in recent yearsETO. The 
first step for studying these systems, has been to employ the spherical jellium model (JM) 
along with the density functional formalism (DFF)QilJ. However, the spherical- JM despite 
its initial successes suffers from two main deficiencies. 

Firstly, the approximation that metal clusters assume spherical geometry, can only be jus- 
tified for large closed-shell clusters, but not for very small clusters having closed- or open-shell 
electronic configurations. Therefore, some authors have used the deformed spheroidal- or 
ellipsoidal- JMstinla. The results of ab initio molecular dynamics calculations^ confirm the 
overall shapes predicted by the deformed- JMs. Koskinen and coworkers^ in their ultimate- 
JM have assumed the jellium background to be completely deformable both in shape and in 
density. However, the ultimate- JM is applicable only when the local density parameter, r s , 
does not differ much from 4.18, i.e., the revalue for which the jellium system is in mechanical 
equilibrium. In the following we will discuss the mechanical stability of the JM. 

The second drawback arise from the JM itself. It is well-known that the JM yields 
negative surface energies^ at high electron densities (r s < 2), and negative bulk mo duliS 
for r s 6. These drawbacks are expected to manifest themselves in the jellium metal 
clusters too. To overcome these deficiencies, some authors have brought in the ionic structure 



equilibrium, except for r s m 4.18. In particular, we refer to the work of Perdew, Tran, and 
SmithH who have introduced the stabilized jellium model (SJM). Applications of the SJM 
to infinite and semi-infinite simple metalsBEil yield realistic estimates for the cohesive and 
surface energies. In fact, the SJM-calculations of the energetics of simple metal clusters with 
spherical geometries!!, and voids@ show an appreciable improvements over the simple JM 
results. For example, the energy per particle for large sodium metal clusters changes from 
~ —2eV to ~ — 6eV. 

However, according to what was mentioned earlier, the spherical-SJM may only be suit- 
able for clusters with closed-shell electronic configurations. The reasons are two-fold: i) In 
the SJM, the stabilization has been accomplished for a spin-compensated system which is not 
necessarily suitable for open-shell clusters, ii) The electronic charge densities for open-shell 
systems do not have spherical symmetry. 

Thus, to study open-shell clusters, two modifications over the spherical-SJM need be 
done. Firstly, open-shell clusters are not generally expected to have zero polarization. This 
is particularly true for clusters with odd number of valence electrons. Therefore, the SJM 
must be generalized to the stabilized spin-polarized jellium model (SSPJM) which is the 
subject of this work. Secondly, because of the deformations due to the Jahn- Teller effectill 
for open-shell systems, one should employ deformed nonspherical shapes for the jellium 
background. Rothlisberger and Andreonit3 by using the Car-Parrinelld2a method, i.e., the 
unified DFF and molecular dynamics, have obtained interesting results. The results of their 
extensive computer simulation show that on the one hand, the overall shapes of clusters 
change when the number of atoms in clusters change. This supports the deformed jellium 
models. On the other hand, they show that the average distance between the nearest 
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neighbors, d av , in sodium microclusters alternate with increasing number of sodium atoms 
(see Fig. 15 of Ref. [pqQ . We will show that these alternations are predicted via the SSPJM. 

Guided by these results, we have asked whether it is possible, to keep the spherical 
geometry and further improve the results obtained using the spherical-SJM. Our answer 
is positive^. In fact, in the spherical-SSPJM, by allowing the volume of a cluster change 
as a function of its spin-polarization, we arrive at new results which are absent in the 
previous spherical- JM or spherical-SJM and gain further insights for the observed odd- 
even alternations in the ionization potentials (IP) of metal clustersSSS. For clarification 
purposes, we mention that it is the difference of the level separations near the Fermi level for 
the neutral and ionized clusters which gives rise to alternations in IPs. These changes can 
show up: i) By allowing nonspherical shape deformations of the jellium background (JB) 



allowing its volume per valence electron change. The latter is what we are addressing in this 
paper. We will see that the mechanism causing changes in the JB volume is rooted in the 
stabilization of the spin-polarized JM. At any rate, in reality, both effects of nonspherical 
deformations and volume changes of the JB are in operation. 

Our self-consistent SSPJM-calculations of the energetics of metal clusters (Al, Ga, Li, 
Na, K, Cs) show that the total energy of a cluster is minimized for a configuration with 
maximum spin compensation (MSC). That is, for clusters with even number of electrons, 
the energy minimization gives rise to complete compensation (iVj = iVj, and for clusters 
with odd number of electrons, only one electron remains uncompensated (iVj — = 1). 
According to our calculations, the only exceptions to the MSC-rule for both neutral and 
ionized metal clusters (Al, Ga, Li, Na, K, Cs) of various sizes (2 < iV < 42) are AI12, Al^ 3 , 
Al 14 , and Gai2- Here, N is the total number of valence electrons and N/z = n is the number 
of z-valent atoms in a cluster. The MSC-rule together with the monotonically increasing 
variation of the cluster size as a function of spin polarization give rise to the alternation 
of the mean distances between the nearest neighbors as a function of cluster size, N. As a 
result of these alternations, the total energies change, and thereby, the IPs alternate. In this 
work, by taking a diffusion layer for the jellium^, we have also repeated our calculations for 
Na-clusters. Application of the diffuse-stabilized spin-polarized jellium model (dif-SSPJM) 
has lowered our calculated ionization energies, and as a result, we have obtained closer 
agreement with experimental data. 

It would be interesting to compare the results of our calculations, i.e., the MSC-rule, with 
the results of the previous calculations in which the density parameter, r s , was assumed to be 
fixed for all cluster sizes. The situation in those models is similar to that of atoms in which 
the external potential produced by the positive charge background with given revalue is 
fixed for different electronic configurations and only the electrons are allowed to redistribute 
themselves. We note that in the case of atoms, clearly the external potential due to nuclear 
charge remains fixed. In these cases according to the Pauli exclusion principle, electrons 
with parallel spins are kept apart, i.e., further apart as compared with the electrons having 
antiparallel spins. In fact, this is the way the electrons reduce their total electrostatic energies 
and this is why the electrons in a given shell assume maximum polarization consistent 
with the Pauli exclusion principle. In other words, the Hund's first rule is applicable in 
these cases. On the other hand, in the case of the SSPJM, the total energy of the system 
is reduced by allowing the cluster radii to expand and the ions as well as electrons to 



while keeping its volume 




spherical shape of the JB and 
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redistribute themselves. As we mentioned earlier, this results in the applicability of the 
MSC-rule instead of the Hund's first rule. 

The organization of this paper is as follows. In section we formulate the SSP JM , and 
show that in view of the Pauli exclusion principle the equilibrium bulk density parameter, r s , 
for nonzero polarization is somewhat different from that of the spin-compensated system. 
Section [Hi] is devoted to the calculational scheme and the application of the SSPJM to 
metallic clusters. In section [TV], we present the results of our calculations of total ground state 
energy per valence electron, the surface and curvature energies, the ionization potentials, 
and the energy second difference, A 2 (iV), which is a measure of relative stability, for different 
metal clusters. Lastly, we explore the condition of the equality of the Fermi energies for the 
up- and down-spin bands (e F = ep) in Appendix |A|. 



In this section we generalize the SJM to include a uniform electron system with nonzero 
constant spin polarization, (. This formulation, in the limit of ( = reduces to the SJM. 
The development of the formulation of the SSPJM parallels the formulation of the SJM. 
Here, the non-interacting kinetic energy functional, T s , as well as the XC-energy functional, 
E xc , depend on the polarization, (, of the system. But the classical Coulomb interaction of 
the electrons with pseudo-ions depends only on their relative distances and is independent 
of the spin polarization of the electrons. In a homogeneous system with nonzero electron 
polarization, the total electron density is the sum over the two spin density components, 
i.e., n = ri| + n^. By defining polarization as ( = (rif — n^/n, each component is expressed 
in terms of £ and n 



As in the SJM, we use the Ashcroft empty core pseudopotential" for the interaction 
between an ion of charge z and an electron at a relative distance r: 



where the core radius, r c , will be fixed by setting the pressure of the system equal to zero. 
The average energy per valence electron in the bulk, with density n and polarization (, is 



II. THE STABILIZED SPIN-POLARIZED JELLIUM MODEL 



n T = |(l + C)n 
n i = |(1 - On. 



(1) 




(2) 



£(n, C) = t s (n, C) + e xc (n, () + w R (n, r c ) + e M {n) + e h . 



(3) 



where 



t s (n, = i c *[(i + 5/3 + (i - 5/ V /3 



(4) 



e xc (n, = ^^[(1 + C) 4/3 + (1 - 4/ V /3 + e c (n, Q 



(5) 
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All equations throughout this paper are expressed in Rydberg atomic units. Here t s and 
e xc are the kinetic and exchange-correlation energies per particle respectively. For e c we use 
the Perdew-Wang parametrizationB wr is the average value of the repulsive part of the 
pseudopotential (wr = innr^), and Em is the average Madelung energy of a collection of 
point ions embedded in a uniform negative background of density n (em = — 9z/5r ). All 
nonuniformity of the true electron density n(r) is contained in the band-structure energy 
Ebs- ro being the radius of the Wigner-Seitz sphere, is given by ro = z 1 ^ 3 r s . We note that 
for monovalent metals z — 1, and if one sets z* — 1 for polyvalent metals, one obtains 
reasonable agreement with experiment (see Ref.||33||). In the latter case, each Wigner-Seitz 
cell will be replaced by z smaller cells with volume 47rrf/3 per cell. As in the SJM, we 
assume that is negligibly small compared to other terms in Eq. (|^). Since energy, and 
thereby, pressure in this formalism depend on ( as well as r s , the stabilization of the bulk 
spin-polarized system with given r s - and (^-values, forces the pseudopotential to assume a 
core radius appropriate to these values, i.e., r c = r c (r s , (). In order to stabilize a bulk system 
with equilibrium density n, and constant (, we should set the pressure equal to zero: 

This equation fixes the core radius, r c , as a function of n and £. Using 




t s (r s ,() = --t s (r s ,0 (8) 



C 



r.. 



and 



e x (r a ,C) = ~e a .(r„C), (9) 



equation flj]) results in 



( d \ 9 
2t s (f s ,()+e x (f s ,()-f s ^—j E c {f s X) + -^l + E M {f s ) = Q. (10) 

The solution of this equation at equilibrium density, n, reduces to the following equation for 
r c , which will now depend on f s and £, 

r c (r s ,C) = ^y- |-2t s (f s ,C) - £*(f s , C) + r s e c (r s ,C) -£m(^)J • (H) 

Because of the (^-dependence of r c , the difference potential (<fo)ws becomes polarization- 
dependent. Here (<fo)ws is the average of the difference potential over the Wigner-Seitz cell 
and the difference potential, 5v, is defined as the difference between the pseudopotential of 
a lattice of ions and the electrostatic potential of the jellium positive background. As in the 
SJM (see Eq. (27) of Ref. |33[|), at equilibrium density we have through Eq. (0) 
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(Sv) ws = n 



d_ 

dn 



[w R (n) + e M (n)) 



—n 



dn 



{t s (n,()+£xc(n,()) 



(12) 
(13) 



so that 



(5v)w = -n(J^J (t a (n,C) + +e x (n,C)) + j(J^j e c (f s ,() (14) 



= ~ lutein, C) + e x (n, C) - r s ^ J^j £c{r s , C) } ■ 



(15) 



In the above equations, n = 3/(47rff) is the equilibrium electronic density of a homoge- 
neous system which has a nonzero constant polarization. Note that, here, n is polarization- 
dependent. In fact, by increasing the polarization, one increases the number of the spin-up 
relative to the spin-down electrons, and therefore, as a consequence of the Pauli exclusion 
principle, the total number of Fermi holes corresponding to the spin-up electrons is increased. 
This leads to the volume expansion of the system ( It is a well-known fact that in a molecule, 
bond lengths depend on spin configurations). In order to estimate the changes of equilibrium 
density as a function of polarization, we put the pressure of the homogeneous spin-polarized 
free electron gas equal to zero. But for a homogeneous spin-polarized free electron gas with 
a uniform positive background, we have 



e(r a , = t a (r 9 , () + e x (r s , () + e c (r s , (), 



(16) 



where the electrostatic energies have cancelled out. Now, vanishing of the pressure at equi- 
librium leads us to 



= I -t s (r s , C) + ye x (r a , () - (JpJ e c (r s , () } , 



A7vrl | r s 



(17) 



the solution of which gives the equilibrium r s as a function of (. At £ — 0, the latter 
equation yields the well-known paramagnetic value of 4.18, and therefore, in the case of 
non-zero polarization, it is convenient to write the equilibrium r s -value as 



ff G (C) =4.18 + Ar s EU . 



EG 



(18) 



The increment Arf G is used to find a rough estimate for the equilibrium revalues of various 
simple metals through 



rf(C) = rf(O) + Arf- , 



EG 



(19) 



where the superscript X refers to a given simple metal. In applying this result to different 
simple metals, we assume that this increment is independent of the value of f s (( = 0), and 
simply depends on the value of (. This is the simplest assumption we have thought of. Other 
forms of ff(() is possible. For instance, one could use a low degree polynomial of ( with 
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coefficients depending on the type of metal. For Al, Ga, Li, Na, K, and Cs, f s (0) are taken 
to be 2.07, 2.19, 3.28, 3.99, 4.96, 5.63 respectively. The increment of equilibrium density n 
due to increasing (, affects the value of the core radius of the pseudopotential and causes 
it to increase monotonically as a function of (. Also the value of (5v )ws will depend on 
(. Fig. [I] shows the behavior of (<5u)ws as a function of polarization, (, for the six metallic 
systems considered in here. One notes that for metals with the value of f s (0) less than 4.18, 
the increase in ( weakens the effective potential relative to its value at ( = 0, but for metals 
with f s (0) > 4.18, the depth of the potential increases. Once the values of (<5v)ws an d r c as 
a function of f s and ( are found, the equation (23) of Ref. can be generalized to 



+~ / dr 4>([n,n + ];r)[n(r) -ra+(r)], (21) 



E S spM n h n h n +] = E m [n h n h n + ] + (e M (n) +w R (n,()) J dr n+(r) 

+ (Sv) ws (n,C) J <2re(r)[n(r)-n + (r)], (20) 

where 

Em[n h ni,n + } = T s [n T ,nJ + E xc [n h n;] 
1 
2 

and 

ttln,n4;r) = J dl ' W\-_ n ;W . (22) 

O(r) has the value of unity inside the system and zero outside. By taking the variational 
derivative of -Esspjm with respect to spin densities nf,nj,, one finds the Kohn-Sham (KS) 
effective potentials 

< f f(rX) = j^j(E SSPm -T s ) (23) 

= 0([n,n + ];r) +< c (r) + e(r)(^) ws (n, C), (24) 

where a =|, \- The distinction between the v"jj for the two schemes of the SSPJM and the 
SJM is rooted in the ^-dependence of the quantities n + and (<5f)ws for the SSPJM. The 
forms of v° c (r) are the same in both cases. By solving the KS-equations 

(V 2 + < // (r))0-(r)=^0-(r) ; a=ti (25) 
n(r) = Yl n *(*)> ( 26 ) 

<T=U 

n„{v) = £ | | 2 , (27) 

i(occ) 

and finding the self-consistent values for el and 0°", one obtains the total energy. In the next 
section we apply this model to simple metal clusters. 
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III. CALCULATIONAL SCHEME 



In the spherical-SSPJM with sharp boundary, the positive background charge density 
is constant and equals to n(C) inside the sphere of radius -R(C) = ^ 1//3 ^s(C) an d vanishes 
outside. In the case of dif-SSPJM, since there exists no sharp boundary for the jellium, we 
have taken the effective boundary at R(() = A^ 1 / 3 f s (<^). Here, for metal clusters, we take 
the value of ( as 

( = (N 1 -N l )/N ; N = N 1 + N l , (28) 

where iVj and N± are the total numbers of valence electrons with spins up and down respec- 
tively. One could instead take £ as the following average: 



C = - / Anr 2 dr y lw ^ n ; = — r , (29) 
Jo n(rj 3 

in which ri is the radius at which the density components fall to the value of, say 1 percent 
of their peak values. In the homogeneous case, the equation for ( reduces to our assumed 
relation. The densities in the integrand of ( are the self-consistent values obtained by solving 
the KS-equations. It turns out that these values of £ and ( are very close to each other and 
so we use the simpler one. The effective potentials in the KS-equation will be 

<//(r; [?., CD = v b (r; [f J) + v n (r; [f s ]) + < c (r; [f s , C]) + e(r)(^) ws (f s , C), (30) 

where in the case of sphere with sharp boundary, 

, v j -(N/R)[3-(r/R) 2 ] ; r<R 
< r) = \ -2N/r -r>R. (31) 

In the above equation, Vj, is the potential energy of interaction between an electron and the 
positive charge background which depends on ( via r s (£). Also we have 

v u (v) = 2 / J^Lr dv', (32) 
J | r — r | 

_ 6E xc [n hni ] 

V * c{T) - 6n a (r) ■ [66) 
A spherical jellium with a finite surface thickness (diffuse-jellium)^l, is defined by 

( \ - i H 1 ~ ( R + t)e- R ^[smh(r/t)]/r} ; r < R . ■ 

n+[T) ~\n{l- ((R + t)/2R)(l -e- 2R /*)}Re( R - r V t /r ; r > R, [ ] 

where R = A^ 1/,3 r s , and t is a parameter related to the surface thickness (for other forms of 
diffuse-jellium see Ref.||). In our numerical calculations with diffuse-jellium for Na-clusters, 
we have chosen t = 1.0 both for neutral and singly ionized clusters. In the case of dif-SSPJM 
the potential energy of an electron due to the background will be 

dif ( \ _ f (47r/3)n(3i? 2 - r 2 - 6t 2 ) - 8nnt 2 (R + t)e- R ^[smh(r/t)]/r ; r<R . . 
Vb [T) ~ \ -2N/r + Annt 2 [(R + t) + (R- ^e^e-^^/r ; r > R. [ } 

Here, R is the effective radius of the jellium sphere, i.e., R = N 1 / 3 ?^). 



8 



IV. RESULTS AND DISCUSSIONS 



In this section we discuss the calculated results for metal clusters (Al, Ga, Li, Na, K, 
Cs) of different sizes (2 < N < 42). After computing the KS-orbitals for spin-up and 
spin-down components, and finding the corresponding eigenvalues, the total energy of an 
A-electron cluster, E(N), were calculated for both sharp- and diffuse-jellium spheres. For 
the sake of comparison, we have also repeated the calculations based on the JM and the 
SJM using the local spin density approximation (LSDA). For a given A-electron cluster, 
we have considered the various combinations of and iVj_ values while keeping the total 
number of valence electrons, N = + N^, fixed. 

Our calculations based on Eq. (0), show that the total energy of the system decreases as 
its polarization decreases, i. e., the system moves towards spin-compensated configurations. 
In fact, the total energy minimization is accomplished when electron spin compensation is 
maximum. Namely, for clusters with even number of electrons, the energy minimization gives 
rise to complete compensation (iVj — N^), and for clusters with odd number of electrons, 
only one electron remains uncompensated (Aj- — = 1). 

In order to show the effect of MSC-rule, we have plotted AMSC= -E^pped — -Emsc, as 
a function of N in Fig. || for Na- and Al-clusters. Here, -Emsc is the total energy of a 
cluster assuming MSC-configuration, and -Empped is the total energy of the same cluster, but 
with a configuration involving only one spin-flip in the outer shells relative to the MSC- 
configuration. The spin-flip is to be consistent with the Pauli exclusion principle. In fact, 
we have calculated AMSC for Al n , Ga n , Li n , Na n , K n , and Cs n -clusters with 2 < nz < 42, 
and found that it is always positive except for Ali 2 , A\± 3 , Ali 4 , and Gai2- Thus, we conclude 
that MSC-rule is at work for nearly all clusters considered in here. 

It is worth noting that in the JM and the SJM calculational schemes the minimum energy 
corresponds to configurations with maximum spin polarization in the outermost shells; i.e., 
Hund's first rule is at work. For example, in a 13-valence electron cluster, the JM and the 
SJM result in the minimum energy configuration in which all 5 electrons in the outermost 
shell (I = 2) are in parallel spin- up state (iVj — N± = 5), while in the SSPJM the energy is 
minimized when — iVj_ = 1. 

The reasons for different behaviors resulting from the use of the JM and the SJM or 
of the SSPJM and the dif-SSPJM, namely the applicability of the Hund's first rule or the 
MSC-rule lies in the unnecessary constraint of rigid jellium background assumed in earlier 
JMs. Note that, in these models one fixes the r s - value for all cluster sizes with arbitrary spin 
configurations. Lifting such a contraint is consistent with the well-known fact that molecular 
bond-lengths depend on spin polarization. Moreover, from molecular dynamics calculations 
one can infer the alternating variation of the nearest neighbor distances as a function of 
cluster sizeJH, which supports our idea of allowing the alternating volume expansion of the 
jellium background positive charge distribution. In cases where the outermost shell is closed 
(N = 2, 8, 18, 20, 34, 40, . . .), all the four schemes of the JM, SJM, SSPJM, and dif-SSPJM 
predict the same spin configurations and the radius of the jellium sphere is the same for 
the first three schemes, but differs in the dif-SSPJM. Also, the total energy values are the 
same both in the SSPJM and the SJM. In cases where the outermost shell contains only 
one electron or lacks one electron to have a closed shell (N = 1, 3, 7, 9, 17, 19, 21, 33, 35, 
39, 41, . . .), the above-mentioned four schemes predict the same spin configurations, but the 
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energies are all different. Here, the difference in the SJM- and the SSPJM-values arise from 
different jellium radii in the two schemes for a given cluster size. In the above two special 
cases, the Hund's first rule and the MSC-rule are identical. 

When the cluster size remains fixed while the polarization changes, as was assumed in 
the case of the JM and the SJM, the only way the system could reduce its total energy was 
to redistribute its electrons further apart. One finds the situation in the JM and the SJM 
to be very similar to that of atoms in which the external potential of the nucleus is kept 
fixed for different spin configurations. As in the case of atoms, due to the Pauli exclusion 
principle, parallel spin electrons are required to stay apart, i. e., further apart than when 
they assume a spin-compensated configuration. In the SSPJM considered in this paper, the 
relative positive charge background radii, R(N, ()/R(N, 0), of clusters are allowed to expand 
with increasing polarization. In a sense, here, the ionic motions are simulated through such 
an expansion. Now, because of the above-mentioned freedom, as soon as electrons try to 
take advantage of the Pauli exclusion principle and begin to spill out of the cluster, the 
positive charge background will try to follow them. In other words, the freedom of cluster 
size expansion renders the application of Hund's first rule unnecessary, and the SSPJM 
chooses to be in a spin-compensated configuration. In short, contrary to the case of the 
spherical- JM and the spherical-SJM, which are governed by the Hund's first rule, here the 
MSC-rule is in effect, and it is this MSC effect which results in the changes in the level 
separations near the Fermi level for neutral and ionized clusters, and thereby, gives rise to 
the well-known odd-even alternations in IPs of alkali metal clusters. One notes that our 
viewpoints about the changes in level separations, and those that attribute this effect to 
nonspherical shape deformations refer to two complementary effects. That is, in our case, 
the shape has remained spherical whereas the volume is allowed to change, but in the case of 
deformed- JMs, the volume is fixed and the shape is allowed to change. In a combined effect, 
one allows, at the same time, the volume of a cluster to change and its shape to deform. 

The size of jellium sphere in the SSPJM for neutral and singly ionized clusters are differ- 
ent. This should be contrasted with the cases of the JM and the SJM in which the sizes for 
the neutral and ionized clusters are assumed to be the same. Our self-consistent calculations 
show that the polarization (, for metal clusters in their minimum energy configuration as 
a function of total number of valence electrons, N, satisfies (with exceptions for Al 12 , Alf 3 , 
AI14, Gai2) the following equation: 

C(N) _ [1 ~ (-1)*] _ f ;iVeven 

U7Vj _ 2N ~ \ 1/N ; N odd. W 

Using the above equation for polarization, we plot the equilibrium revalue for different 
Na-clusters in Fig. |||. As seen from the figure, for an Na-cluster with even number of atoms, 
r s equals 3.99 and with odd number of atoms the envelope is a decreasing function of N. 
According to Eq. ([19|) , Fig. |3| for other metals remains the same but shifted according to 
the value of r s (0). As is expected, in the limit of large- iV clusters, addition or removal of 
an electron does not change the configurations of all other ions. Comparing Fig. ^ with 
Fig. 15(a) of Ref. [KSJ for average nearest neighbor distances, one notes the same staggering 



effect. 

Fig. |]a and Fig. |]b show the values of the nonbulk binding energies, E/N — ay, for 
Na- and Al-clusters using the JM, the SJM, and the SSPJM along with the LSDA. Here, 
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ay represents the bulk binding energy. The calculated energies based on the JM, the SJM, 
and the SSPJM for Na-clusters, aside from the details, are nearly close and positive as 
expected. This is because the sodium revalue (3.99) is close to the zero-pressure jellium 
r s -value (4.18). However, in the case of Al-clusters, the energies based on the SJM and 
the SSPJM remain near each other and are positive, but the energies based on the JM is 
considerably away from them and at times assume negative values (see Fig. f|b). The reason 
is that the JM at high densities fails (note that the r s -value of Al is 2.07) and leads to 
mechanical instability. We have also calculated the nonbulk binding energies of other metal 
clusters (Ga, Li, K, Cs) and see that as r s decreases, the energies based on the JM moves 
away from the results based on the SJM and the SSPJM. By looking at Fig. [1] we note that 
the (<k>)ws contribution to the total energies are positive in the case of Cs and K (for which 
r s > 4.18); and negative in the case of Al, Ga, Li, and Na (for which r s < 4.18). Thus, in 
the case of K- and Cs-clusters we should expect that the JM energies using the LSDA lie 
above the SJM- and the SSPJM-energies, and in the case of Na-, Li-, Ga-, and Al-clusters 
we should expect it to lie below them. Our calculations of the nonbulk energies of A1-, Ga-, 
Li-, Na-, K-, Cs-clusters confirm these conclusions. 

In Fig. | we show the plot of A 2 (A r ) = E(N + 1) + E(N - 1) - 2E(N) which determines 
the relative stability of different Na-clusters. We have also calculated A 2 (iV) for other metal 
clusters mentiond in this paper. In the plot using the SSPJM, the incorrect peaks predicted 
by the spherical-SJM at iV =5, 13, 27, 37 have disappeared and the clusters with N =8, 18, 
20, 40 are predicted to be more stable. These results are consistent with the fine-structure 
observed experimantally in the abundance curved. Similar observations can be made from 
the plots of A 2 (iV) for other metals considered in here. The overall agreement between our 
results and experimental data is good. 

Next, we have calculated the surface and curvature energies of A1-, Ga-, Li-, Na-, K-, 
Cs-clusters. According to the liquid drop mo one may write the total energy of a 

finite quantal system in the form of 

E = E + 5E (37) 

where 5E is shell correction and E is the smooth part of the total energy which is written 
as a sum of volume, surface, and curvature contributions. In the case of spherical geometry, 
E reduces to the following parametrized equation as a function of the number of valence 
electrons in a neutral cluster 

E(N) = a v N + a s N 2/3 + a c N 1/3 . (38) 

Here, ay is the total energy per electron in the bulk. Its absolute values for the metals 
Al, Ga, Li, Na, K, Cs using the SJM or the SSPJM are respectively 10.57, 10.14, 7.37, 
6.26, 5.18, 4.64 in units of electron-volts using the r s - values mentioned earlier. The surface 
energy, a, and the curvature energy, Ac, are related to the parameters as and ac through 
the relations: 

1 A 1 

a s , A c = a c . (39) 



The parameters as and ac are obtained by a least-square fit of our self-consistent total 
energies to Eq. (|38f) . 
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Fig. ^a shows the surface energies of different metals as a function of their revalues for 
the three schemes using the SSPJM, the SJM, and the JM, and compared with the results 
obtained by others (BULK)0. The surface energies of the SSPJM and the SJM remain 
positive for the metals considered here, but it becomes negative in the case of the JM (as 
mentioned in the introduction) for high electron density metals (Al, Ga). 

Fig. |6|b shows the curvature energies of various metals as a function of their revalues. 
The SSPJM results in a higher curvature energies than the other two schemes using the SJM 
and the JM. The results of bulk calculations^ are also shown. 

Finally, in Figs. [7]a-f, we have compared the ionization energies of metal clusters for 
different schemes as well as with the experimental values. The ionization energy of an N- 
valence electron cluster is defined as the difference between the total ground state energy of 
that system with N and (N — 1) valence electrons. The experimental values are taken from 



Ref. |49| for aluminum, Ref. |4fl for lithium, sodium, and potassium. When we calculate the 
IPs by means of the SSPJM and the SJM, we see that the odd-even alternations in alkali 
metal clusters, not present in the results of the JM or the SJM results, show themselves up 
in the SSPJM calculations. In the SJM using the LSDA, kinks appear only at half-filled 
and closed shells, whereas in both the SSPJM and the dif-SSPJM, there exists a peak at 
the middle of each pair of adjacent odd numbers. These peaks correspond to clusters with 
even number of atoms, in agreement with experimental data. In dif-SSPJM, the ionization 
energies of Na-clusters are lowered so that the values near the closed shells (where the 
nonsphericity becomes less important) are in good agreement with experiment than the 
corresponding values of the other two models. At the end, because of the assumed spherical 
geometry, the pronounced shell effects are still present in the IPs when we go from one closed 
shell cluster to a cluster with one more electron and the overall saw-toothed behavior of IPs 
remains. In other metals the agreement between theory and experiment is poor for all the 
three schemes. We think our results for Na-clusters will improve if dif-SSPJM is used along 
with spheroidal or ellipsoidal geometries. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have generalized the SJM to the spin-polarized case by allowing the 
volume of the spherical positive charge background to change and have calculated the en- 
ergetics of metal clusters (Al, Ga, Li, Na, K, Cs). Our self-consistent calculations show 
that for spherical geometries, the minimum-energy of a metal cluster is obtained when the 
electronic spin compensation is maximum. That is, in contrast to the spherical- JM and 
the spherical-SJM which are governed by the Hund's first rule, here the MSC-rule is in 
effect. We have discussed in detail that the situations in the JM and the SJM are similar 
to that in atoms. In both cases, the external potential - being due to the positive charge 
background or due to the nuclear charge - are fixed and therefore, in both cases the Hund's 
first rule is applicable. However, in the case of the SSPJM, because of the extra degrees 
of freedom, namely the expansion of the positive charge background, the system assumes 
maximum spin-compensated configurations. This MSC-rule results in the alternations of 
the average distance between the nearest neighbors, and thereby, in the alternations of the 
IPs. Moreover, application of the dif-SSPJM for alkali metal clusters brings the IP-values 
closer to the experimental data. Finally, we believe that if the SSPJM is used in conjunction 
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with nonspherical shape deformations, better agreements between theory and experiment 
will result. 
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APPENDIX A: STABILITY CONDITIONS IN A SPIN-POLARIZED INFINITE 

JELLIUM 

In an unpolarized infinite jellium (£ = 0), the stability of the system is obtained via 
de/dr s \^=o— 0. In other words, in this case, it suffices to establish mechanical stability. 
This is done by employing the Ashcroft empty-core pseudopotential whose core radius is 
fixed by setting the pressure equal to zero. However, for a spin-polarized infinite jellium, 
the total energy per electron, e, depends on n and (. Thus 

& (i)„ d <' < A1 > 

At equilibrium, the following two equations must be satisfied: 

= (A2) 




Eq. ( |A~2| ) establishes the mechanical stability and Eq. ( |A3| ) sets the Fermi energy of the up- 
and down-spin bands equal, i. e., e f = s F . The latter point is proved below. e\ is the 
highest occupied spin-up KS-level and ep is the highest occupied spin-down KS-level. 
proof: 

Following the work of Russier, Salahub, and Mijouleil, we can write 



de\ de ( dn^\ de ( dn^ 



\d( J dn^ \ d( J n dni \ d( 
so that using Eq. (JJ) reduces to 



(A4) 



(|). = (A5) 

Now, using the variational principle for the ground state energy of a homogeneous system 
subject to the two constrains 

J n a dr = N a , a =|, |, (A6) 
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we have 



S {J^ + n £ ( n h n i) dv ~ M j n T dr - VlJ n l dr ] = °- ( A7 ) 



This equation results in the following Euler equations: 



de 

e + n- fxt = 0, (A8) 



de , , , 



Subtracting Eq. ( |A"9f) from Eq. (|A8|) and dividing both sides by n, one obtains 



y dn-i dni J n 
Substituting Eq. ( |A10| ) into Eq. ([A5|) , one obtains 

(de\ 1 



de de \ 1 . . . . > 

-(M-Vi)- (A10) 



, A ; 2 (MT — A*i)- (AH) 



Eqs. ( [All ) and ( |A3[ ) in conjunction with Koopmans' theoremS entails the equality of the 



two Fermi levels for an infinite homogeneous spin-polarized jellium. 

However, in this paper we have calculated the total energy of finite clusters. In such 
systems, because of the discrete nature of the energy eigenvalues, the equality of the spin- 
up and spin-down Fermi energies does not hold any more. Thus, the stabilization will be 
established through Eq. (|A2| ) and when dE/d( \ Ts changes sign. In practice, for a given (, 
we have first calculated the core radius entering the Ashcroft pseudopotential via Eq. (|TTD 
and used it as an input for the total energy calculation of the cluster. We then, varied (, 
i. e., the cluster's spin configuration, till total energy minimization was attained. Although, 
the procedure employed here is a consistent choice, but one could start with a two parameter 
pseudopotential and apply Eq. (|A2|) and Eq. (|A3|) to the stabilized spin-polarized infinite 
jellium. In this way, one obtains the dependence of the two mentioned parameters on ( 
and r s and repeats the clusters' total energy calculations again. Work in this direction is in 
progress. 
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FIGURES 



FIG. 1. Equilibrium potential difference in Rydbergs versus polarization for Al, Ga, Li, Na, 
K, and Cs metals. 

FIG. 2. The behavior of AMSC as a function of number of valence electrons, N, in units of eV. 
Both neutral and singly ionized metal clusters (Al and Na) are presented. AMSC is the difference 
between the energy of the maximum spin compensated electronic configuration from the energy of 
a different configuration in which one electron is flipped. In the case of Na-clusters all AMSCs are 
positive and in the case of Al-clusters, except for AI12, Al^, AI14, all AMSCs are again positive. 

FIG. 3. Equilibrium Wigner-Seitz radius in atomic units as a function of number of atoms in 
an Na-cluster. 

FIG. 4. The nonbulk total energies per atom of (a) Na- and (b) Al-clusters as a function of 
number of valence electrons in the cluster. The symbols, solid squares, diamonds, and triangles 
correspond to the SSPJM, the SJM, and the JM respectively. For Al, the physical points (i.e., 
multiples of 3) are specified and for Na, the large square symbol corresponds to dif-SSPJM. 

FIG. 5. The second difference of total energies of Na-clusters in units of eV as a function of 
number of atoms in a cluster. The symbols, solid squares, triangles, and diamonds correspond to 
SJM, JM, and dif-SSPJM respectively. 

FIG. 6. (a) Surface energies of metals, a, as a function of their revalues, in units of erg/cm 2 ; 
(b) Curvature energies of metals, A c , as a function of their revalues, in units of mili-hartree/bohr. 
The symbols, crosses, solid squares, triangles, and upside-down triangles correspond to the SSPJM, 
SJM, JM, and bulk calculations (BULK) (see Ref. |Q). The negative surface energies of Al and 
Ga are due to the instability of the JM. 

FIG. 7. The ionization energies in units of eV as a function of the number of atoms in a 
cluster for (a) Al, (b) Ga, (c) Li, (d) Na, (e) K, (f) Cs. The large solid squares represent the 
experimental data. The small solid squares, the diamonds, the triangles, and the upside-down 
triangles, respectively correspond to calculational schemes based on the SSPJM, the SJM, the JM, 
and the dif-SSPJM. By moving towards closed shell electronic structures, the agreement between 
the results of the dif-SSPJM and experimental data improves. 
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